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The IceCube Neutrino Observatory has observed a diffuse flux of TeV-PeV astrophysical neutrinos 
at 5.7cr significance from an all-flavor search. The direct detection of tau neutrinos in this flux has 
yet to occur. Tau neutrinos become distinguishable from other flavors in IceCube at energies above 
a few hundred TeV, when the cascade from the tau neutrino charged current interaction becomes 
resolvable from the cascade from the tau lepton decay. This paper presents results from the first 
dedicated search for tau neutrinos with energies between 214 TeV and 72 PeV in the full IceCube 
detector. The analysis searches for IceCube optical sensors that observe two separate pulses in a 
single event - one from the tau neutrino interaction, and a second from the tau decay. No candidate 
events were observed in three years of IceCube data. For the first time, a differential upper limit on 
astrophysical tau neutrinos is derived around the PeV energy region, which is nearly three orders 
of magnitude lower in energy than previous limits from dedicated tau neutrino searches. 
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I. INTRODUCTION 

The IceCube Neutrino observatory has announced a 
significant detection of a diffuse astrophysical neutrino 
flux above 30 TeV [U [2] . The source of this flux is as yet 
unknown, with possible candidates including cosmic ray 
acceleration in active galactic nuclei, gamma ray bursts, 
and supernova remnants. Assuming standard 3-flavor os¬ 
cillations and a most commonly considered \ \ 

= 1 : 2 : 0 flux from pion decay at the source, the 
neutrinos detected in IceCube should be divided almost 
equally into all three flavors [3]. Other flavor composi¬ 
tions at the source ranging from 1 : 0 : 0 to 0 : 1 : 0 
are possible for dominant processes such as neutron de¬ 
cay a, energy loss of pious and muons before decay in 
environment with strong magnetic fields or high matter 
density IMH] and muon acceleration [9]. Though those 
scenarios result in non-universal flavor ratios at Earth, 
they all predict significant fluxes of tau neutrinos after 
averaged oscillations by propagation over astronomical 
distances [iDHia. Above PeV energies, the Earth be¬ 
comes opaque to electron and muon neutrinos, while the 
tau neutrino flux is regenerated through subsequent tau 
lepton decays to neutrinos [13] . Tau neutrino background 
from the atmosphere is expected to be negligible at high 
energies, with only a small contribution from the decay 
of charmed mesons [T4|. Therefore the detection of tau 
neutrinos at high energies would both give new infor¬ 
mation about the astrophysical flux as well as serving 
as an additional confirmation of the astrophysical ori¬ 
gin of the high energy diffuse neutrino signal. Two re¬ 
cent flavor ratio analyses of IceCube high energy neutrino 
events were consistent with equal fractions of all flavors 
in IceCube, though with large uncertainty dHHl]. Nei¬ 
ther flavor ratio analysis included a dedicated tau neu¬ 
trino identification algorithm, which would improve the 
measurement of astrophysical neutrino flavor ratios. Pre¬ 
cise measurement of astrophysical neutrino flavor content 
at Earth will shed light on the emission mechanisms at 
the source, test the fundamental properties of neutrinos 
over extremely long baselines and better constrain new 
physics models which predict significant deviations from 
equal fractions of all flavors [T5]-[24]. 

Most neutrino interactions in IceCube have one of two 
event topologies: tracks from charged current (CC) in¬ 
teractions of muon neutrinos, and cascades (or showers) 
from CC interactions of electron and low energy tau neu¬ 
trinos and neutral current (NC) interactions of all fla¬ 
vors. As the average tau decay length roughly scales as 
5 cm/TeV, at energies above a few hundred TeV, the tau 
lepton produced in a tau neutrino CC interaction would 
have a decay length sufficiently long that the CC interac¬ 
tion of the tau neutrino and the subsequent decay of the 
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FIG. 1. Sketch of a Vr undergoing a charged current interac¬ 
tion and producing a double bang signature in IceCube. The 
blue circles represent IceCube photon sensors buried in the 
ice. This figure is not drawn to scale. 

tau lepton could be resolved by IceCube sensors. There 
is an 83% chance that the tau lepton decays to electrons 
or hadrons, producing a second cascade. This double cas¬ 
cade signal is called a “double bang” |3] ; an event sketch 
is shown in Eigure[^ 

A previous IceCube search for high energy tau neu¬ 
trinos was published using data from the partially com¬ 
pleted detector [25]. Since a fully contained double bang 
is unlikely to be contained in the small size of the in¬ 
complete detector, the search was optimized for partially 
contained double bangs. The search led to a null result 
and was in fact more sensitive to electron and muon fla¬ 
vor neutrinos than to tau neutrinos. 

There is a 17% chance of the tau decaying to a muon, 
producing an outgoing track. Therefore, another possi¬ 
ble search method is to look for a track that abruptly 
brightens along its length as the muon produces light 
more efficiently than the parent tau [26]. Tau neutrinos 
may also produce a signature in cosmic ray air shower 
detectors. The Pierre Auger Observatory has reported 
an upper limit on the tau neutrino flux from a search 
for horizontal showers from tau lepton decays induced 
by Earth-skimming cosmogenic neutrinos. However, the 
energy threshold of this search is 200 PeV, much higher 
than the energy of the astrophysical neutrinos observed 
by IceCube [27] . 

This paper describes a new search method for double 
bangs whose two cascades may not be separately recon¬ 
structed, but which appear as a two-peaked or “double 
pulse” waveform in one or more IceCube sensors. 

II. THE ICECUBE DETECTOR 

The IceCube detector [28] consists of 86 vertical cables, 
called strings, deployed in the ice near the geographic 
South Pole. Each string contains 60 Digital Optical Mod¬ 
ules (DOMs). A DOM consists of a 10 inch photomulti¬ 
plier tube (PMT), digitizing electronics, and LED flash¬ 
ers for calibration [29]. The digitized PMT signal is called 
a waveform. The DOM utilizes two digitizers: the Ana¬ 
log Transient Waveform Digitizer (ATWD) which digi- 
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tizes at 3.3 ns per sample for 128 samples, and the fast 
Analog to Digital Converter (fADC) which digitizes con¬ 
tinuously at 25 ns per sample, with a record length set 
at 256 samples. The ATWD output is separated into 3 
different gain channels (xl6, x2, xO.25) to cover the dy¬ 
namic range of the PMT, which has a linear response 
(within 10%) up to currents of 400 photoelectrons (PE) 
per 15 ns m- When a PMT receives a signal above a 
threshold of 0.25 PE, this is called a hit. The xl6 gain 
channel is captured first, with the x2 and xO.25 chan¬ 
nels captured if the next lowest gain channel exceeds 768 
ADC counts in any sample. A local coincidence hit (LC) 
occurs if a pair of nearest or next-to-nearest neighbor 
DOMs on the same string are hit within 1 microsecond. 
Eor LC hits, the complete ATWD and fADC waveforms 
will be sent to the surface. The primary IceCube trigger 
keeps all DOM hits if 8 or more LC hits occur anywhere 
in the detector within a 5 microsecond window; such a 
collection of DOM hits is called an event. 

IceCube employs a number of filtering algorithms in 
order to reduce the data volume for transmission to the 
Northern hemisphere. The analysis described in this pa¬ 
per uses the “Extremely High Energy” (EHE) filtering 
algorithm, which keeps all events that deposit more than 
1000 PE in the detector. 

IceCube was fully built as of December 2010. This 
analysis uses 914.1 days of data from the full detector 
between May 13, 2011 and May 6, 2014. The data were 
kept in the analysis chain only when all IceCube strings 
were operating and no in-situ calibration light sources 
were in use. 

Background and signal passing rates were computed 
using Monte Carlo simulation. The CORSIKA [ST] sim¬ 
ulation package is used to generate cosmic ray induced 
muons. Astrophysical and atmospheric neutrinos are 
simulated using an adapted version of the Monte Carlo 
generator ANIS [32]. Photon propagation through the 
ice is simulated as described in [33]. PMT response and 
digitization electronics are fully simulated, which is par¬ 
ticularly important for this analysis. 

In addition to the simulation, 10% of the data were 
used to develop cuts and estimate cosmic ray muon back¬ 
ground rates, with the rest of the data not used until the 
cuts were finalized. 


III. SEARCH FOR TAU NEUTRINOS 
A. Double Pulse Algorithm 

The goal of the double pulse algorithm (DPA) is to 
identify double pulse waveforms that are consistent with 
Ut CC interaction signatures in IceCube, while reject¬ 
ing waveforms with features that are consistent with late 
scattered photons from single cascade events from NC 
and Ue CC interactions. There are two additional types 
of background events which could produce substantial 
double pulse waveforms: (1) high energy single muons 



Time [ns] 

FIG. 2. A simulated double pulse waveform obtaiued iu oue 
DOM from a Ur CC eveut. The primary ueutriuo euergy for 
this eveut is 2.4 PeV, about 75% of this euergy trausfers to 
the outgoiug t~ leptou, which travels 40 meters before decay. 
Iu this eveut, a total of 34 double pulse waveforms were pro¬ 
duced from adjaceut DOMs ou ueighboriug striugs uear the 
eveut iuteractiou vertices. The distauces from the CC vertex 
aud the t~ decay vertex to the DOM that produced this dou¬ 
ble pulse waveform are 76 m aud 75 m, respectively. Time 
= 0 correspouds to the begiuuiug of the eveut readout wiu- 
dow, which begius 10 microsecouds before the eveut trigger 
lauuches. 

and/or muon bundles induced by cosmic rays interacting 
with the atmosphere; (2) CC interactions in IceCube 
which produce energetic muons. For double pulse wave¬ 
forms caused by energetic atmospheric muons, the first 
pulse is usually from a combination of Cherenkov light 
emissions and coincident stochastic energy loss, and the 
second pulse is from TeV-scale stochastic energy losses 
tens of meters away from the DOM. For double pulse 
waveforms from astrophysical CC events, the first 
pulse is from energy deposition of the CC hadronic inter¬ 
action vertex, while the second pulse is from a coincident 
stochastic energy loss of the energetic outgoing muon. 
Since the double pulse waveforms from CC events, 
energetic atmospheric muons and astrophysical CC 
events are not distinguishable from one another as they 
are caused by the same mechanism of two substantial en¬ 
ergy depositions near certain DOMs, we do not remove 
these events with the DPA. They are to be removed later 
by comparing their overall topologies and timing profiles. 
The potential impact of instrumental backgrounds on the 
DPA was found to be negligible. Afterpulses in the PMT 
waveforms, caused by ionization of residual gases by elec¬ 
trons accelerated in the space between dynodes, usually 
occur from 500 ns to microseconds later than the primary 
pulse, and therefore do not appear as double pulses in a 
single waveform. Late pulses, caused by photoelectrons 
backscattered from the first dynode, occur on a time scale 
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of 60 ns later than the primary pulse, but usually have a 
low amplitude and do not trigger the DPA [29] . 

The DPA uses the positive and negative first deriva¬ 
tives of a waveform to determine rising and trailing edges. 
A double pulse is defined as a rising edge, followed by 
a trailing edge, followed by another rising edge. Wave¬ 
forms from the ATWD digitizer in the lowest gain chan¬ 
nel available are used since higher gain channels are usu¬ 
ally saturated for high amplitude waveforms. The fADC 
waveforms are not used since they do not have multi¬ 
ple gain channels available and since their coarser timing 
causes double pulse features to be blended together or 
saturated. The DPA uses 7 configurable parameters to 
characterize a double pulse waveform: 

• Since signal waveforms are from bright events close 
to a DOM, the DPA is only run on ATWD wave¬ 
forms that have integrated charge greater than qi 
= 432 PE. 

• The beginning of the waveform is determined by 
a sliding time window of 3.3 ns which searches for 
a monotonic increase in the waveform amplitude 
within a time span of 3.3 x 6 = 19.8 ns. 

• Once the beginning of the waveform is identified, 
the waveform is divided into segments of 4 ATWD 
bins (13.2 ns) and the first derivative of the wave¬ 
form is computed in each segment. 

• If the first derivative is positive in ni = 2 consecu¬ 
tive segments, this is considered the rising edge of 
the first pulse. When the subsequent derivative is 
negative for 77-2 = 2 consecutive segments, this is 
considered the trailing edge of the first pulse. The 
rising edge of the first pulse is required to have an 
integrated charge of at least q 2 = 23 PE, and the 
integrated charge of the trailing edge is required to 
be at least qs = 39 PE. The integrated charge sums 
up all the charge corresponding to the entire rising 
or trailing edges, which usually last longer than two 
segments (26.4 ns) for a large pulse. 

• The second pulse rising edge is defined when the 
derivative after the trailing edge of the first pulse 
is positive again for 77,3 = 3 consecutive segments. 
This requirement is due to the fact that the light in 
the second pulse is often more scattered and there¬ 
fore has a less steep rising edge than the first pulse. 
The second pulse trailing edge is often outside the 
ATWD window, and hence is not included in the 
calculation. The rising edge of the second pulse is 
required to have an integrated charge of at least q^ 
= 42 PE. 

The configurable DPA parameters rii and qi were 
tested and optimized using a variety of IceCube event 
waveforms including simulated neutrinos of all flavors, 
simulated atmospheric muons and data from in-situ laser 
calibration devices. An example of a simulated double 



FIG. 3. Fraction of events that pass the charge cut of 
logio(QTot)>3.3 and have at least one double pulse wave¬ 
form as a function of total deposited charge. The higher the 
charge in the event, the more likely it is to contain at least one 
double pulse waveform. Muons with lower deposited charge 
that might evade containment cuts are less likely to produce 
double pulse waveforms. 


pulse waveform from a CC event is shown in Fig- 
ure[^ The production of double pulse waveforms depends 
largely on the distances between event interaction ver¬ 
tices and nearby DOMs. The median distance between 
the Ur CC (tan lepton decay) vertices and the double 
pulse DOMs is 49 (44) meters. 

The fraction of events that pass a charge cut of 
logic (QTot) >3.3 and have at least one double pulse wave¬ 
form is shown in Figure as a function of deposited 
charge, for both the 10% data sample and atmospheric 
muon simulation. Near the charge cut threshold, fewer 
than 1 in 1000 events will include a double pulse wave¬ 
form. 


B. Event Selection 

The event selection process was carried out at three cut 
levels, driven by the specific goal of background rejection 
at each level. To conform with standard IceCube usage, 
the cuts for this analysis are numbered beginning with 
Level 4. The cut levels are summarized as follows: 

Level 4: Events are required to have at least one wave¬ 
form which passes the DPA. An additional event-wise 
charge cut of logic(QTot)>3.3 is also required to enrich 
the sample with high energy events. 

Level 5: At this stage, we remove track-like double 
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FIG. 4. Events per year that pass the level 6 containment cut on depth and distance to the detector edge (solid line) for 
simulated atmospheric muons (left) and astrophysical Vr (right). 


pulse events which are predominantly due to atmospheric 
muons. Following the Level 4 cut, each event is recon¬ 
structed using a maximum likelihood method based on 
a hypothesis of an infinite track and a hypothesis of a 
point-like cascade. These reconstructions only make use 
of the timing information for the earliest photon arriving 
at the DOMs and hence are computationally efficient. 
The reduced log likelihood ratio between the two hy¬ 
potheses Li?=log(Lcascade/Ltrack) IS required to be neg¬ 
ative, indicating the event topology is more cascade-like 
than track-like. This cut eliminates most down-going en¬ 
ergetic muons and muon bundles. To further veto down¬ 
going muons, the first hit in the event is required to be 
below the top 40 meters of the instrumented volume. 
CORSIKA simulation predicts that 3.5T3.4 atmospheric 
muons survive to Level 5 in 914.1 days. 

Level 6: At this stage we eliminate cosmic ray in¬ 
duced muons which pass near the edges of the detector 
and hence appear cascade-like. An additional reconstruc¬ 
tion algorithm is performed on all events which pass the 
preceding cuts, using full charge and time information, 
which is more computationally expensive. A boundary is 
defined by the surface connecting the position of the out¬ 
ermost layer of strings in the detector. The containment 
criterion requires that the reconstructed vertex be inside 
the instrumented volume and a given distance away from 
the boundary. The distance from the boundary depends 
on depth, with stricter containment required at the top 
and bottom of the detector. The containment is illus¬ 
trated in Figure]^ which shows the distribution of event 
vertices with respect to the boundary for signal and for 
atmospheric muon background. Due to the scarcity of 
events in the atmospheric muon simulation at high en¬ 
ergies, the double pulse criterion is removed from this 
plot, with all other cuts kept. The very few atmospheric 
muon events which survive the containment cut are close 
to the charge cut threshold, and as shown in Figure 
such events have a lower probability of producing a dou- 



FIG. 5. Passing rate for signal (astrophysical Vr GG events in 
solid blue) and backgrounds (astrophysical in dot-dashed 
magenta, astrophysical Ve in dot-dashed green, atmospheric 
neutrinos in dotted purple, and atmospheric muons or GOR- 
SIKA in dashed red) as a function of cut level. Data shown 
here in solid black is 10% of the total data sample. 


ble pulse waveform. 

At the final cut level, the predicted rates from all 
sources in three years of data are summarized in Ta¬ 
ble [ij The assumed astrophysical flux is based on the 
diffuse flux measured by IceCube at the level of 
= 1.0 X 10“^ GeV cm“^ s~^ sr“^ per flavor [2], and 
90% of the predicted CC events are between 214 TeV 
and 72 PeV. A softer neutrino energy spectrum of E“^-^ 
[12] reduces the expected number of Ur events and the 
dominant background of CC events by 36% and 57% 
respectively. The atmospheric neutrino rate prediction 
includes both neutrinos from n/K decay [34| and neutri- 
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TABLE L Predicted event rates from all sources at the final 
cut level. Errors are statistical only. 


Data samples 

Events in 914.1 days (hnal cut) 

Astrophysical Ur CC 
Astrophysical CC 

Astrophysical Ue 
Atmospheric ly 
Atmospheric muons 

(5.4 ± 0.1) • 10“^ 

(1.8 ± 0.1) • 10“^ 

(6.0 ± 1.7) • 10“2 
(3.2 ± 1.4) • 10“^ 

(7.5 ± 5.8) • 10“2 


nos from charmed meson decay mi. The primary cosmic 
ray spectrum used to predict atmospheric neutrino rates 
is corrected for air shower measurements in the knee re¬ 
gion of several PeV [35] . 

Figure summarizes the passing rate of signal and 
background events at each cut level. At the final cut 
level, astrophysical events have the highest passing 
rate of any source, and the dominant background is as¬ 
trophysical CC events. The effective areas for CC 
and CC events at the final cut level is shown in Fig¬ 
ure An optimal energy window for the astrophysical 
Vr search in IceCube using this double pulse method is 
around the PeV region, where the effective areas for 
CC events are nearly an order of magnitude higher than 
that of CC events. It is planned that events found at 
final cut level will be further investigated with segmented 
energy loss reconstruction algorithms [36] to acquire their 
energy loss profile and directionality. Event probabilities 
of z^T-like or not will also be computed based on likelihood 
methods. 



FIG. 6. Effective areas at final cut level as a function of pri¬ 
mary neutrino energy. Only the middle 90% Vt energy range 
(214 TeV - 72 PeV) is plotted. The dominant background 
for this analysis is due to astrophysical CC events, so only 
Vt CC (solid blue) and CC (dashed red) effective areas 
are shown. The plot demonstrates that the optimal energy 
window for the astrophysical Vr search using the double pulse 
waveform approach is from 0(100) TeV to 0(10) PeV. In par¬ 
ticular, around PeV energies, effective areas for Vr CC events 
are about an order of magnitude higher than those for CC 
events. Effective areas for Ve (not shown) are 1-2 orders of 
magnitude below the effective areas for CC, except at the 
Glashow resonance energy of 6.3 PeV m- 


IV. SYSTEMATIC UNCERTAINTIES 

The sources of systematic uncertainty considered in 
this analysis are neutrino cross sections, anisotropy in 
the optical scattering in ice, uncertainty in the optical 
scattering and absorption lengths in ice, and DOM ef¬ 
ficiency. The main sources of systematic uncertainty in 
the signal are summarized in Table [ll] 

The neutrino cross sections used in this analysis are 
from the CTEQ5 model [38]. The CSMS model [39] . 
which has updated parton distribution functions, pre¬ 
dicts ^ 5% fewer events compared to the CTEQ5 model. 

An earlier study in IceCube attempting to reconstruct 
the double deposition of energy from a z^^- CC event has 
found that the recently identified anisotropy in the op¬ 
tical scattering in ice m would modify the number of 
expected photons in some DOMs and hence could mimic 
a double cascade feature in the reconstructed energy seg¬ 
ments [41] . A study based on simulations with and with¬ 
out this anisotropy found a 7% lower signal event rate 
prediction for the double pulse analysis when anisotropy 
was included. The effect is small at the waveform level 
due to the fact that the double pulse events usually oc¬ 
cur within tens of meters of a DOM, which is within 1-2 


scattering lengths in the ice. 

The optical scattering length and absorption length 
were varied according to the uncertainty in the value of 
these parameters [33|. Increasing the absorption by the 
allowed uncertainty decreases the signal event rate by 
4.9%, and decreasing the absorption and scattering in¬ 
creases the signal event rate by 8.1%. 

Since the Ur double pulse events are very bright, uncer¬ 
tainty in the DOM efficiency does not play an important 
role. Simulation with the DOM efficiency set at -1-10% 
and -10% of the nominal values yielded a decrease of 
1.6% in the signal event rate when decreasing the effi¬ 
ciency, and an increase of 6.7% in the signal event rate 
when increasing the efficiency. 

Adding the various errors in quadrature, the total sys¬ 
tematic uncertainty in the signal is about ±10%. 

The uncertainty in the atmospheric muon and neu¬ 
trino background is dominated by statistical error, due 
to the fact that few simulated background events pass 
the cuts. The largest source of systematic error is un¬ 
certainty in the cosmic ray flux at high energies which 
contributes ±30%/-50% uncertainty to the atmospheric 
muon flux and ±30% uncertainty to the atmospheric neu¬ 
trino flux m- 






















FIG. 7. Event 1 before level 6 containment cut with its corresponding double pulse waveform. This event occurred on May 
30, 2011. The colored spheres indicate hit DOMs, with size indicating the amount of charge deposited on the sphere and color 
indicating time: red is earlier, blue is later. 




FIG. 8. Event 2 before level 6 containment cut with its corresponding double pulse waveform. This event occurred on November 
27, 2011. 


f 



FIG. 9. Event 3 before level 6 containment cut with its corresponding double pulse waveform. This event occurred on August 
28, 2012. 
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TABLE IL Source of systematic uncertainty in the signal. 


Neutrino cross sections 

-5% 

Anisotropy in the optical scattering in ice 

-7% 

Optical scattering and absorption lengths in ice 

+ 8 . 1 % 

- 4 . 9 % 

DOM efficiency 

+ 6 . 7 % 

- 1 . 6 % 

Total 

+ 10 . 5 % 
- 10 . 0 % 


RESULTS 


Zero events were found after all cuts were applied. At 
level 5, before the containment cut, three events were 
found which each have one double pulse waveform, all of 
which occurred on strings at the edge of IceCube. These 
events are consistent with atmospheric muons interacting 
near the edge of the detector, producing a double pulse 
waveform in a cascade-like event but failing the subse¬ 
quent containment cut at Level 6. The observation of 3 
events in 914.1 days of livetime matches the CORSIKA 
prediction at level 5 as discussed in Section \nrB\ The 
events and their corresponding double pulse waveforms 
are shown in Figures and 

Based on zero observed events, an integrated astro- 
physical Ur flux upper limit is set to be = 5.1 x 

10“^ GeV cm“^ sr“^ s“^. A Ur flux differential upper 
limit in the energy range of 214 TeV to 72 PeV, which 
contains 90% of the predicted CC events, is also ex¬ 
tracted following the procedure that was employed in de¬ 
riving quasi-differential upper limits from previous EHE 
cosmogenic neutrino searches in IceCube [milHlIlH]. In 
this procedure, flux limits were computed for each en¬ 
ergy decade with a sliding energy window of 0.1 decade, 
assuming a differential neutrino flux proportional to 
jSO]. Since zero events were found, the 90% C.L. 
event count limit in each energy decade is 2.44 based 
on the Feldman-Cousins approach m- The dominant 
sources of systematic error in this analysis are indepen¬ 
dent of energy. Therefore, all the sources of systematic 
and statistical error are incorporated in the limit cal¬ 
culation by uniform scaling of the effective area. The 
differential upper limit is plotted in Figure pT| 


VI. CONCLUSION 


The double pulse search method is shown to be robust, 
with the observed background from cosmic ray induced 
muons matching prediction. The search is more sensi¬ 
tive to tan neutrinos between 214 TeV and 72 PeV than 
to any other flavor. Given the astrophysical neutrino 
flux observed by IceCube, fewer than one tan neutrino 
candidate event is expected in three years of IceCube 
data, and none are observed. A differential upper limit 
has been placed on the astrophysical tan neutrino flux. 



FIG. 10. Neutrino flux upper limits and models as a function 
of the primary neutrino energy. The thick red curve is the Ur 
differential upper limit derived from this analysis, including 
systematic and statistical errors. In computing the differential 
upper limit, values of the flux limit were calculated for each 
energy decade with a sliding energy window of 0.1 decade. 
The thick black error bars depict the all-flavor astrophysical 
neutrino flux observed by IceCube [2] . The thick dashed line is 
the differential upper limit derived from a search for extremely 
high energy events which has found the hrst two PeV cascade 
events in IceCube |42l |43] . The blue dotted line is the Auger 
differential upper limit from Ur induced air showers [27] . The 
orange dashed line is the Waxman-Bahcall upper bound which 
uses the UHECR flux to set a bound on astrophysical neutrino 
production EH. The dash-dotted line (magenta) represents 
the prompt neutrino flux predicted from GRBs; prompt in 
this context means in time with the gamma rays [45]. The 
dash-dot-dot line (grey) indicates the neutrino flux predicted 
from the cores of active galaxies (JG). The thin dash-triple¬ 
dot line (red) shows the neutrino flux predicted from starburst 
galaxies, which are rich in super novae m- 


with an energy threshold three orders of magnitude lower 
than previous dedicated tan neutrino searches by cos¬ 
mic ray air shower detectors. Searches for double bang 
events with well separated cascades in IceCube are under¬ 
way. Future extensions of IceCube such as the proposed 
IceCube-Gen2 detector [52] will have a factor of 5 to 10 
times more sensitivity to astrophysical tan neutrinos than 
the current IceCube detector. 
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